Strain-engineered magnetic order in (LaMnOa) n / (SrMn03)2„ superlattices 
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Using first-principles calculations based on the density functional theory, we show a strong strain 
dependence of magnetic order in (LaMn03) n /(SrMn03)2n (001) superlattices with n — 1,2. The 
epitaxial strain lifts the degeneracy of Mn e g orbitals, thus inducing an inherent orbital order, which 
in turn strongly affects the ferromagnetic double exchange of itinerant e g electrons, competing with 
the antiferromagnetic superexchange of localized t2 g electrons. For the case of tensile strain induced 
by SrTiOa (001) substrate, we find that the ground state is A-type antiferromagnetic and d x 2_ y 2 
orbital ordered, which is in excellent agreement with recent experiments [S. J. May et ai, Nature 
Materials 8, 892 (2009)]. Instead, for the case of compressive strain induced by LaAlC>3 (001) 
substrate, we predict that the ground state is C-type antiferromagnetic and d 3z 2_ r 2 orbital ordered. 

PACS numbers: 75.70. Cn, 71. 20. Be, 75.47.Lx, 75.25.Dk 



I. INTRODUCTION 

Transition metal oxides in perovskite based structures 
exhibit a wide variety of phases with different electronic, 
magnetic, and orbital structures, and show rich function- 
alities such as high-Xc superconductivity, colossal magne- 
toresistance, and multiferroicsi A recent advance in epi- 
taxial growth techniques has made it even possible to fab- 
ricate transition metal oxide heterostructures with sharp 
and smooth interfaces controlled at the atomic scaled 
In these heterostructures, many unique properties, not 
found in the corresponding alloy compounds made of the 
same composite elements, have been observed, which in- 
clude e.g., two dimensional electron gas with high mo- 
bility at the heterostructure interfaces^ indicating the 
promising potential of oxide heterostructures for future 
technological applications^ 

In the case of manganites^ LaMnOs is an A-type an- 
tiferromagnetic insulator and SrMnC>3 is a G-type an- 
tiferromagnetic insulator. On one hand, the randomly- 
cation-doped alloy Lai-xSr^MnOs exhibits a rich mag- 
netic phase diagram, depending on the doping con- 
centration x. On the other hand, La/Sr cation- 
ordered analogues forming superlattices behave quite 
differently from their alloy compounds^— For exam- 
ple, La2/3Sr!/ 3 Mn03 alloy has a mixed valence of 
Mn 3+ /Mn 4+ , and the ground state is ferromagnetic half 
metallic due to the double exchange mechanism^ To the 
contrary, it is found experimentally that cation-ordered 
(LaMn0s)2n/(SrMn03) n (001) superlattices are insulat- 
ing when n is larger than 3i&£ This change of behavior 
is easily understood because the number n of SrMn03 
layers control the quantum confinement potential: when 
n is small, the confinement potential is small and the e g 
electrons are distributed uniformly, thus expecting the 
phases similar to the alloy Lai-^Sr^MnOs. When n is 



large, the confinement potential becomes large enough to 
trap the e g electrons in LaMnOs layers, and thus the bulk 
properties of LaMnOs and SrMnOs would be observed. 
Several theoretical studies for (LaMn03)2n/(SrMn03)„ 
superlattices have been reported to understand their elec- 
tronic and magnetic properties ! 11 ^ 12 

More recently, Bhattacharya et a/. 13 ' 14 have experi- 
mentally studied the transport and the magnetic prop- 
erties of similar superlattices (LaMn03)„/(SrMn03)2 n 
grown on SrTi0 3 (001) substrate. They have found that 
the ground state of these superlattices with n = 1,2 
are A-type antiferromagnetic metals with Neel temper- 
ature (Tn) which is higher than that observed in any 
alloy Lai-zSr^MnO;} compound^ Although the simi- 
lar physical principles found in (LaMn03)2n/(SrMn0 3 )„ 
superlattices are certainly expected to apply here, 
the systematic theoretical investigations are required 
to understand the main ingredients which determine 
the electronic as well as the magnetic properties of 
(LaMnOs)™/ (SrMnOs)2n superlattices. 

Here, in this paper, performing first-principles cal- 
culations based on the density functional theory, we 
study the electronic and the magnetic structures of 
(LaMn0 3 )„/(SrMn0 3 ) 2n (001) superlattices with n = 
1,2. We show that the magnetic properties are governed 
not only by the quantum confinement potential caused 
by periodic alignment of cation ions La 3+ /Sr 2+ but also 
by the strain induced by substrates on which the super- 
lattices are grown. Namely, for the case of tensile strain 
induced by SrTi0 3 (STO) (001) substrate, our calcula- 
tions show that the ground state of these superlattices 
are A-type antiferromagnetic and d x i_ y i orbital ordered 
with higher Tn for n = 1 than for n = 2. This is indeed 
in excellent agreement with recent experimental obser- 
vations^ Instead, for the case of compressive strain in- 
duced by LaA103 (LAO) (001) substrate, we predict C- 
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type antifcrromagnctic and d 3z 2_ r 2 orbital orders with 
higher Tn for n = 1 than for n = 2. 

The rest of this paper is organized as follows. After de- 
scribing the computational details in SecHH the numeri- 
cal results for the cases of SrTiC>3 substrate and LaA10 3 
substrate are presented in Sec. MI Al and Sec. IIII B[ re- 
spectively, followed by discussion of the confinement po- 
tential in Sec. IIII CI Sec. II VI summaries this paper. 

II. COMPUTATIONAL METHODS 

We perform the first-principles electronic structure cal- 
culations based on the projected augmented wave pseu- 
dopotcntials using the Vienna Ab initio Simulation Pack- 
age {VASP)^& The valence states include 3p4s3d and 
2s2p for Mn and O, respectively. The electron interac- 
tions are described using the generalized gradient approx- 
imation (GGA) and the rotationally invariant GGA+ U 
method^— with the effective U c g, i.e., U — J, from 1 
cV to 5 cV for d electron states. Compared to the GGA, 
the GGA+J7 approach gives an improved description of 
d electron localization^ The atomic positions of super- 
lattices are fully optimized iteratively until the Hellman- 
Feynman forces are 0.01 eV/A or less. The plane- wave 
cutoff is set to be 500 eV and al2xl2xl2 Monkhorst- 
Pack k-point grid is used in combination with the tetra- 
hedron method^ 

The supercells considered here consist of 6 Mn02 lay- 
ers, 2 LaO layers, and 4 SrO layers for both n = 1 
and 2, as shown in Fig. Q] (a) and Fig. [5] (a). We con- 
sider 12 and 10 different magnetic moment alignments 
to search for the ground state magnetic structures for 
LaMn0 3 /(SrMn0 3 ) 2 and (LaMn0 3 )2/(SrMn0 3 )4 super- 
lattices, as shown in Fig. [3] and Fig. [H respectively. 
These magnetic structures include not only simple ferro- 
magnetic, A-type, C-type, and G-type antiferromagnetic 
structures^ but also magnetic structures with mixed 
combinations of these simple magnetic structures. The 
epitaxial constraint on these supcrlatticcs, which is grown 
on substrates, is to fix the in-plane lattice constants. 
Thus, to simulate the strain effect, we fix the in-plain 
lattice constants (a) of the superlattices to the ones of 
substrates, i.e., a = 3.905 A for SrTi0 3 substrate^ and 
a = 3.81 A for LaA10 3 substrate^! and the lattice con- 
stant (c) perpendicular to Mn02 layers is fully relaxed. 
Atomic positions are also fully optimized. 

III. RESULTS 

A. (LaMn0 3 )„/(SrMn0 3 ) 2 „ on SrTiO :i 

Let us first examine (LaMn0 3 ) n /(SrMn0 3 ) 2 „ (001) 
superlattices on SrTi0 3 (001) substrate. Our system- 
atic GGA calculations reveal that the ground states of 
these superlattices with n = 1 and 2 are both A-type 
antifcrromagnctic metals. A schematic spin alignment 



of A-type antifcrromagnctic order is shown in Fig. [3] (g) 
and Fig. 2] (a). Indeed, as shown in Fig. [1] (d) and Fig. [2] 
(d), the projected spin density distribution, calculated 
by integrating spin density of occupied states from Fermi 
level down to —0.5 eV, clearly indicates the A-type anti- 
ferromagnetic spin order. Our GGA+U calculations also 
find that these A-type antifcrromagnctic states arc ro- 
bust against electron correlations, and they are indeed 
stable up to U c s — 2 eV for n = 1 and U e e = 1.3 eV for 
n = 2 (see Fig. E])^ 

Since the supercell sizes and the numbers of each type 
of atoms are the same, we can simply compare the total 
energy of these two different superlattices. Tab. U sum- 
marizes the total energies for the A-type antiferromag- 
netic states and other magnetic states. Since the A-type 
(C-type) magnetic structure is ferromagnetic (antiferro- 
magnetic) within the ab plane and antiferromagnetic (fer- 
romagnetic) along the c direction, we can approximately 
estimate an effective magnetic exchange ( J e ff) simply by 
comparing the total energy of the A-type and the C-type 
antifcrromagnctic states. It is clearly observed in Tab. U 
that the stabilization energy of the A-type antiferromag- 
netic state, i.e., J e ff, is larger for n = 1 than for n = 2. 
This implies that Tn for n = 1 is higher than that for 
n = 2. These results are in excellent agreement with 
experimental observations by May et al^ 

TABLE I. Total energies (in unit of eV) of 
(LaMn03) n /(SrMn03)2n superlattices (n = 1,2) cal- 
culated using GGA. FM, A-AFM, and C-AFM stand 
for ferromagnatic, A-type antiferromagnetic, and C-type 
antiferromagnetic states, respectively. 



n 


SrTiOs substrate 


LaA103 substrate 


FM 


A-AFM 


C-AFM 


FM 


A-AFM 


C-AFM 




-466.262 


-466.636 


-465.854 


-465.852 


-465.459 


-466.466 


n=2 


-465.637 


-465.885 


-465.151 


-465.398 


-465.015 


-465.858 



Since the epitaxial constraint of substrates is to fix 
the in-plane lattice constant a of the superlattices, the 
tetragonal distortion should inevitably occur, which in 
turn affects the relative occupation of Mn e g electrons. 
Indeed, as shown in Tab. UH we find that the SrTi0 3 
substrate induces tensile strain with a > c, in which 
d x 2_ y 2 orbital is lower in energy than d 3z 2_ r 2 orbital. 
This can be seen in the projected charge density distri- 
bution, the integrated charge density from Fermi level 
down to -0.5 eV, shown in Fig. [T] (b) and Fig. [2] (b), 
indicating that e g electrons preferably occupy d x 2_ y 2 or- 
bital. Because of this orbital order induced inherently by 
the substrate strain, the A-type antiferromagnetic order 
is stabilized. Remember that the magnetic interaction 
between Mn ions is determined by competition between 
the ferromagnetic double exchange via itinerant Mn e g 
electrons and the antiferromagnetic superexchange be- 
tween localized Mn ti g electrons. When d x 2_ y 2 orbital 
is occupied rather than d 3z 2_ r 2 orbital, the strong dou- 
ble exchange induces ferromagnetic order in the ab plane 
while the weak itineracy of d x 2_ y 2 electrons along the c 
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FIG. 1. (Color online), (a) A schematic figure of the supercell considered for LaMn03/(SrMn03)2 (001) superlattices, and 
the projected charge [(b) and (c)] and spin [(d) and (e)] density distributions (integrated from Fermi level down to —0.5 eV 
using GGA) for LaMn0 3 /(SrMn0 3 )2 superlattices grown on SrTi03 [(b) and (d)] and LaAlOs [(c) and (e)] (001) substrates. 
The loci of MnO layers are indicated in (a), where red, blue, green, and purple spheres indicate O, La, Sr, and Mn atoms, 
respectively. In (d) and (e), the up and down spin densities are denoted by yellow and light blue, respectively. 




FIG. 2. (Color online), (a) A schematic figure of the supercell considered for (LaMnOs)2/(SrMn03)4 (001) superlattices, and 
the projected charge [(b) and (c)] and spin [(d) and (e)] density distributions (integrated from Fermi level down to —0.5 eV 
using GGA) for (LaMn0 3 )2/(SrMn0 3 )4 superlattices grown on SrTi0 3 [(b) and (d)] and LaA10 3 [(c) and (e)] (001) substrates. 
The loci of MnO layers are indicated in (a), where red, blue, green, and purple spheres indicate O, La, Sr, and Mn atoms, 
respectively. In (d) and (e), the up and down spin densities are denoted by yellow and light blue, respectively. 



direction reduces substantially the double exchange and 
as a results the superexchange between t2 g electrons sta- 
bilizes antiferromagnetic order along this direction. Fi- 
nally, it is also interesting to note that the optimized 
lattice constant c for n = 1 is shorter than that for n = 2 



(see Tab. [II]), which is also qualitatively in good agree- 
ment with experimental observations ^ 
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FIG. 3. (Color online). 12 different magnetic structures considered for LaMn03/(SrMn03)2 superlattices: G-AFM (a), C-AFM 
(b), Ml-AFM (c), FM (d), M2-AFM (e), D-AFM (f), A-AFM (g), M3-AFM (h), M4-AFM (i), M5-AFM (j), M6-AFM (k), and 
Dl-AFM (1). Mn spins are indicated by arrows. Aqua, lime, and violet spheres stand for Sr, La, and Mn atoms, respectively. 
O atoms are omitted for clarity. 



TABLE II. The optimized lattice constant c (averaged 
value within the supercell and in unit of A) and c/a of 
(LaMn03)„/(SrMnC>3)2n superlattices (n = 1,2) calculated 
using GGA. The magnetic structures are A-type and C-type 
antiferromagnetic for SrTi03 and LaA103 substrates, respec- 
tively. 



n 


SrTi03 substrate 


LaA103 substrate 


c 


c/a 


c 


c/a 


n=l 


3.806 


0.9746 


4.006 


1.0115 


n=2 


3.825 


0.9795 


4.010 


1.0525 



B. (LaMn0 3 )n/(SrMn0 3 )2n on LaAlOs 

Now, let us study the electronic and the magnetic prop- 
erties of (LaMn03)„/(SrMn03)2n superlattices (n = 
1, 2) on (001) LaAlC>3 substrate. In the alloy manganitcs 
Lai-zSr^MnOs, it is known that c/a is a key parameter 
in determining the magnetic ground states^ Here, we 
demonstrate that even in these superlattices the mag- 
netic structure can be controlled by the substrate strain 
which varies c/a. 

Because the in-plane lattice constant of LaAlC>3 is 
much smaller than that of LaMnOs (bulk lattice param- 
eter is 3. 935 A), it is expected that the LaAlC>3 substrate 
induces compressive strain. In fact, we find in Tab. |TT] 
that the lattice constant c in the superlattices is larger 



than the in-plain lattice constant a. As a result of this 
tetragonal distortion, Mn e g orbitals are split and d 3z 2_ r 2 
orbital is lower in energy than d x 2_ y 2 orbital, which thus 
induces d 3z 2_ r 2 orbital order. A signature of this orbital 
order can be seen in the projected charge density distri- 
butions shown in Fig.[T](c) and Fig.[2](c). Because of this 
orbital order, the magnetic ground state is expected to be 
C-type antiferromagnetic. Considering 10 — 12 different 
candidates for possible magnetic structures as shown in 
Figs. [3] and [31 our GGA calculations find that the ground 
states for n = 1 and 2 are both C-type antiferromagnetic 
metals [Fig. [3] (b) and Fig. 01(b)] . This magnetic align- 
ment can be indeed clearly seen in the projected spin 
density distribution as shown in Fig. [1] (e) and Fig. [2] 
(c). It is also interesting to note that the lattice distor- 
tion along the c direction is less pronounced for the case 
of LaAlC>3 substrate as compared to the case of SrTi03 
substrate. As shown in Fig. [6J Mn-O-Mn angles between 
the nearest layers along the c direction for the superlat- 
tices on LaAlOs substrate is almost 180°, which certainly 
favors the ferromagnetic double exchange along this di- 
rection. We also find that the C-type magnetic structure 
is robust against electron correlations in Mn d orbitals 
up to £/off = 4 eV for n = 1 and U e s = 1-5 eV for n = 2 
(see Fig.H)^ 

As in the case of SrTiOs substrate, we can discuss the 
Neel temperature TV for the C-type antiferromagnetic 
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FIG. 4. (Color online). 10 different magnetic structures considered for (LaMn03)2/(SrMn03)4 superlattices: A-AFM (a), 
C-AFM (b), D-AFM (c), FM (d), M2-AFM (e), G-AFM (f), Ml-AFM (g), M3-AFM (h), M4-AFM (i), and M5-AFM (j). 
Mn spins are indicated by arrows. Aqua, lime, and violet spheres stand for Sr, La, and Mn atoms, respectively. O atoms are 
omitted for clarity. 
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FIG. 5. (Color online). U D a dependence of the relative ener- 
gies (calculated using GGA+U) for various magnetic struc- 
tures (see Fig. [3] and Fig. |4]) compared to A-type antiferro- 
magnetic state for (LaMn03)„/(SrMn03)2n with n — 1 (left) 
and n = 2 (right) on SrTiOs substrate. 



FIG. 6. (Color online). Mn-O-Mn angles between the nearest 
layers along the c direction for the relaxed crystal structures 
(calculated using GGA) for (LaMn0 3 )n/(SrMn0 3 ) 2 n (001) 
superlattices grown on different substrates indicated in the 
figure. The layer positions in the horizontal axis are indicated 
in Fig. [1] (a) and Fig. [2] (a). 



order by calculating the total energy, and the results are 
summarized in Tab. HI Simply by comparing the total 
energies of the C-type and the A-type antiferromagnetic 



states, the difference of which gives a rough estimate of 
an effective magnetic exchange J c g , we find that the sta- 
bilization energy of the C-type antiferromagnetic state, 
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(LaMn0 3 )/(SrMn0 3 )j 
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(LaMn0 3 ) 2 /(SrMn0 3 ) 4 
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be very useful in estimating the confinement potential 
for more complex superlattices in which first-principles 
electronic structure calculations are comnnta.tionallv ex- 



FIG. 7. (Color online). U e g dependence of the relative ener- 
gies (calculated using GGA+J7) for various magnetic struc- 
tures (see Fig. [3] and Fig. [4} compared to A-type antiferro- 
magnetic state for (LaMn0 3 ) n /(SrMn0 3 )2n with n — 1 (left) 
and n = 2 (right) on LaAlOs substrate. 



i.e., J e ff, is larger for n = 1 than for n = 2. This im- 
plies that Tjv for n = 1 is higher than that for n = 2. 
Since (LaMn03)„/(SrMn03)2n superlattices (71 = 1,2) 
on (001) LaAlC>3 substrate have not been studied experi- 
mentally, these results provide the theoretical prediction 
which should be tested experimentally in the future. 



C. Confinement potential 

Finally, let us briefly discuss why the magnetic and or- 
bital ground states found here are spatially uniform, in 
spite of apparent periodic potential modulation caused 
by different ionic charges, i. e., La 3+ in LaMnOs lay- 
ers, and Sr 2+ in SrMn03. As reported in Ref. [l2|, one 
way to estimate the effective potential modulation is to 
evaluate the oxygen Is core energy level. The results 
for (LaMn03)„/(SrMn03)2n superlattices with n = 1 
and 2 are shown in Fig. HI From these figures, we see 
that (i) the potentials are almost the same for both sub- 
strates, and (ii) as is expected, the confinement potential 
becomes larger with n. The calculated charge density 
shows that e g electrons in LaMn03 layers is ^0.2 (0.1) 
more than that in SrMn03 layers for n = 2 (n = 1). 
This suggests that the thickness is still thin enough not 
to confine e g electrons in LaMn03 layers. However, we 
naturally expect that the bulk properties may recover 
far away from interface when n is increased further and 
a metal-insulator transition should eventually occur. 

It is also found that the confinement potential can be 
more easily estimated simply by calculating Madelung 
potential. As shown in Fig. \§\ Madelung potential can 
indeed semi-qualitativcly reproduce the values estimated 
from oxygen Is core energy level. This finding should 
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FIG. 8. (Color online). The variations of the relative oxygen 
Is core energy (calculated using GGA) in each Mn02 layer 
of (a) LaMn0 3 /(SrMn0 3 )2 and (b) (LaMn0 3 )2/(SrMn0 3 )4 
superlattices. Results for SrTi0 3 and LaA10 3 substrates are 
indicated by black circles and red squares, respectively. The 
layer positions are indicated in Fig. [T] (a) and Fig. [2] (a). 




345 
Layer position 



FIG. 9. (Color online). The variation of Madelung potential 
for Mn ions in each Mn02 layer of (a) LaMn0 3 /(SrMn0 3 )2 
and (b) (LaMn0 3 )2/(SrMn0 3 )4 superlattices. Results for 
SrTi0 3 and LaA10 3 substrates are indicated by black circles 
and red squares, respectively. Here, the ideal crystal struc- 
tures with no distortion, and the ideal Mn valency (indicated 
in the figures) with O 2- , La 3+ , and Sr 2+ are assumed. The 
layer positions are indicated in Fig. [T](a) and Fig. [2] (a). Note 
that the sign convention of Madelung potential used here is 
that electrons prefer to locate in LaMn0 3 layers. 
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IV. SUMMARY 

Using first-principles calculations based on the den- 
sity functional theory, we have studied the effects 
of epitaxial strain on the magnetic ground states in 
(LaMn0 3 ) n /(SrMn0 3 )2 n (001) superlattices with n = 
1,2. Our results clearly demonstrate that as in alloy 
manganitcs, even in superlattices, the epitaxial strain in- 
duced by substrates enforces tetragonal distortion, which 
in turn governs the ground state magnetic structure via 
the inherent orbital ordering. We have found that for 
the tensile strain induced by SrTiC>3 (001) substrate, 
the ground state is A-type antifcrromagnctic metal with 
d x 2_ y 2 orbital order. The approximate estimation of an 
effective magnetic exchange suggests that the Neel tem- 
perature T/v of the A-type antiferromagnetic order is 
higher for n = 1 than that for n = 2. These results are 
in excellent agreement with experimental observations^ 
Furthermore, we have predicted that for the compressive 



strain induced by LaAlC>3 (001) substrate, the ground 
state is C-type antiferromagnetic metal with d 3z 2_ r 2 or- 
bital order with higher Neel temperature Tjy for n = 1 
than that for n = 2. These predictions should be con- 
firmed experimentally in the future. 
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